Exact semi-relativistic model for ionization of atomic hydrogen by electron impact 
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We present a semi-relativistic model for the description of the ionization process of atomic hy- 
drogen by electron impact in the first Born approximation by using the Darwin wave function to 
describe the bound state of atomic hydrogen and the Sommerfeld-Maue wave function to describe 
the ejected electron. This model, accurate to first order in Z/c in the relativistic correction, shows 
that, even at low kinetic energies of the incident electron, spin effects are small but not negligible. 
These effects become noticeable with increasing incident electron energies. All analytical calcula- 
tions are exact and our semi-relativistic results are compared with the results obtained in the non 
relativistic Coulomb Born Approximation both for the coplanar asymmetric and the binary coplanar 
geometries. 
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I. INTRODUCTION 

Relativistic (e, 2e) processes have been reviewed both 
from the experimental and theoretical point of view 
[1]. As one deals with atomic hydrogen, the value of 
the parameter Za is very lower than one, where Z is 
the atomic charge number and a is the fine structure 
constant. Therefore, it is convenient and sufficient to 
use approximate wave functions of a mathematically 
simpler structure than the exact analytical wave func- 
tions needed to describe relativistic (e,2e) processes. 
A numerical approach to an exact description of the 
relativistic ionization of atomic hydrogen by electron 
impact could be carried out but we will focus instead on 
an alternative approach that will give nearly the same 
results as the exact description if the condition Za <C 1 
is satisfied. In (e, 2e) processes, relativistic effects are 
important and all electrons (the incident, scattered 
and ejected) can have very high velocities. One has to 
consider many interactions (to name some, retardation 
interaction, magnetic interaction and spin-dependent 
interaction). For atomic hydrogen, many experimental 
and theoretical contributions have been made [2-3]. 
Some were successful but the theoretical situation for 
all set-ups and kinematics is far from resolved, at least 
analytically. Many calculations have resorted to various 
approximations. For example, plane wave models 
[4-7] are successful in the coplanar binary geometries 
[4] and for fast scattered and ejected electrons. The 
first Born approximation (FBA) has been used to 
describe asymmetric geometries at non relativistic 
energies [8-9]. In this approximation, the incident 
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and scattered electrons are described by plane waves 
whereas the ejected electron is treated as a Coulomb 
wave. Many authors extended this approximation to 
the relativistic domain. Das et al [10-11] employed 
a semi-relativistic Sommerfeld-Maue wave functions 
to describe the ejected electron. Jakubafia- Amundsen 
evaluated the first-order transition matrix element Sfi 
using semi-relativistic Coulomb wave functions times 
a free spinor i.e neglecting the relativistic contraction 
of the bound state and approximating the continuum 
Coulomb state by a relativistic Coulomb wave times a 
free spinor. This model did well in predicting integrated 
cross-sections [12] but yielded a value for the absolute 
triple differential cross section (TDCS) too large. For 
the Coulomb approximation, Jakubafia- Amundsen 
argued that one could not neglect the Coulomb potential 
in the treatment of inner-shell ionization of high-Z 
atoms. Agreement with experiment was encouraging for 
intermediate values of Z . The merits and shortcomings 
of this theory have been analyzed in [14]. Thereafter, a 
fully relativistic version was produced [15] which showed 
that the original physical insight was essentially correct. 

In this contribution, we present a theoretical semi- 
relativistic model, the semi-relativistic Coulomb Born 
Approximation (SRCBA) in a closed and exact form 
for the description of the ionization of atomic hydrogen 
by electron impact that is valid for all geometries. 
Indeed, in the non relativistic Coulomb Born Approx- 
imation (NRCBA), a well known integral occurs [16] 
and is usually denoted by /(A). In this article, we 
show that the main contribution to the unpolarizcd 
triple differential cross section (TDCS) in the SRCBA 
corresponding to the ionization of atomic hydrogen in 
its ground state by electron impact comes from this 
term added to relativistic corrections valid to first order 
in Z/c. These relativistic corrections contain a new 
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integral we have denoted J(A) and in the Appendix, 
we give the formal derivation of this integral. To our 
knowledge, it is the first time that such an integral 
is written down analytically. Needless to say that all 
numerical appropriate tests to check the validity of the 
analytical result we have found have been carried out 
with a very good degree of accuracy. It turns out that 
spin effects can be accounted for even at low kinetic 
energies of the incident electron in the case of the 
Ehrhardt coplanar asymmetric geometry [17] where, for 
a given kinetic energy Tj of the fast incident electron, 
a fast ("scattered ") electron of kinetic energy Tf is 
detected in coincidence with a slow ("ejected") electron 
of kinetic energy Tb- These spin effects as well as the 
relativistic effects become noticeable with increasing 
incident electron kinetic energy. 

The organization of this paper is as follows : in section 
II, we present the semi relativistic formalism of (e, 2e) 
reaction and give a detailed account of the various terms 
that contribute to the unpolarized TDCS, in section III, 
we discuss the results we have obtained and we end by a 
brief conclusion in section IV. The formal derivation of 
the integral J(A) is given in the Appendix. Throughout 
this work, atomic units (a.u) are used (h = m e = e = 1) 
where m e is the electron rest mass. 



II. THE UNPOLARIZED TRIPLE 
DIFFERENTIAL CROSS SECTION 

In this section, we calculate the exact analytical ex- 
pression of the semi relativistic unpolarized TDCS in the 
SRCBA corresponding to the ionization of atomic hydro- 
gen by electron impact. The transition matrix element 
for the direct channel (exchange effects are neglected) is 
given by 



Sfi = -i J dt< i/) Pf (xi)<f>f(x2) | V d | tpp^xi)^^) > 

/+oo r- 
dt I dri^ p/ (i,ri)7 ( ° 1) r/'p 4 (*,ri) 
-OO J 

x < <j> f (x 2 ) | V d | 0,(2:2) > 
In Eq. (1), Vd is the direct interaction potential : 



(i) 



v„ = -L-i 



(2) 



ri are the coordinates of the incident and scattered 
electron, r 2 are the atomic electron coordinates, r\ 2 — 
I ri — r 2 I and n =| ri |. The wave function tjj pi (xi) = 
^ p (t,ri) = u(p, s) exp(—ip.x)/^/2EV is the electron 
wave function described by a free Dirac spinor normalized 
to the volume V and <frij(x 2 ) = <fiij(t,r 2 ) are the semi 
relativistic wave functions of the hydrogen atom where 
the index i stands for the initial state, namely the ground 



state and the index / stands for the final state. The quan- 
tity p.x = p^x^ is the Lorentz scalar product. The semi 
relativistic wave function of the hydrogen atom used is 
the Darwin wave function for bound states [18] 



where 



4>i(t,T 2 ) = cxp(-ie b t)(^ (±) (r 2 ) 



^(r 2 ) = (l4-^a.V (2) )«( ± ),A,(r 2 ) 



(3) 



(4) 



is a quasi relativistic bound state wave function ac- 
curate to first order in Z/c in the relativistic correc- 
tions (and normalized to the same order) with ipo be- 
ing the non relativistic bound state hydrogenic func- 
tion. The spinors are such that = (1,0,0,0) T 
and = (0, 1,0, 0) T and represent the basic four- 

component spinors for a particle at rest with spin-up and 
spin-down, respectively. For the spin up, we have 



¥>< + >(r 2 ) = N D 



1 


± cos 6» 2 



=e- r2 (5) 



V T c sin(0 2 ) exp(i</> 2 ) J 
and for the spin down, we have 



p<->(r 2 ) 



N 



D 



\ 



where 



^ sin(0 2 ) exp(-i^ 2 ) 
V -£cos(fl 2 ) J 



N D = 2c/V4^TT 



\/7T 



(7) 



is a normalization constant lower but very close to one. 
The wave function </>/(£, r 2 ) in Eq. (1) is the Sommerfeld- 
Maue wave function for continuum states [18] also accu- 
rate to the order Z/c in the relativistic corrections. We 
have 4>f(t,r 2 ) = cxp(— iEbt)tp p ^ (r 2 ) and 

i'p^fa) = exp(7rr? s /2)r(l + W7B)exp(ipB.r 2 ) 
x {l 4 - 7^ a -^ (2)} lFii-iTjB,!, -i{p B r 2 + p B .r 2 )) 



u{pb, Sb) 
V^EbT 



(8) 



normalized to the volume V. The Sommerfeld parameter 
is given by 



Eb 

Vb = - — 

c z pb 



(9) 
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In Eqs. (3) and (8), a is related to the 7 Dirac matrices We give the final compact form of the Sommerfed- 
[19] and in the standard representation reads Maue wave function 

•=(:;) <«» 

with a =(<r x , <J y ,<J z ) and the matrices a x ,a y , a z are the 
usual Pauli matrices. The matrix differential operator 
a.V is given by 



1 d z d x - id y 

a.V = 



d x + idy —d z 

d z d x - idy 

\d x + idy -dz 



(11) 



ipp^fa) = cxp(7r?7B/2)r(l + ir] B )cvp(ip B .r 2 )^ iF\{-ir) B , 1, -i(pBr 2 + Ps-r 2 ) + (a.p B + p B a.r 2 ) 

XiFi(-if? B + l,2,-i(p B r 2 + p B .r 2 ))}^g^ (12) 

In Eq. (12), the operator a.p B acts on the free spinor u(j>b, s b ) and the operator a.r 2 acts on the spinor part of the 
Darwin function. The direct transition matrix element in Eq. (1) becomes 

c ■ I a "(P/> s /) u(pb,s b ) o f .. . i , ~ 

% = "V Tll^ 7 ^ V2^b^ 7 W ( ifi(^,l,t(p B r + p B .r)l 4 - — (a.p B + PB a.r) 

x iF^irjB + l,2,i(p B r + p B .r)) jp^r) exp(-ip B .r)[exp(iA.r) - 1] 

x ^S{E s +E B -E t -e b )^0eM^ B /2)T{l- lVB ) (13) 
This transition matrix element contains three terms, one of which is given by 

^ = "V ^^SW^^mvM 1 l>^r + PB.r)))^ >(r) 

87T^ lid)' S '1 

x cxp(-z Pi 3.r)[cxp(2A.r) - 1]^(£/ + £ s - - e t ) exp(^ B /2)r(l - zrys) (14) 

This term can be recast in the form : 

S,< - - i [i/ l(q _A-p B )-H 1 (q-- PB )]^_ 7(1)71 _ 7(2)7 _ — 

x S(E f + E B -S J -e b )exp(7rr/ B /2)r(l-ir ?i j) (15) 
In the above expression, -ffi(q) is given by 

#i(q) = J drcxp(ici.r)iFi(ir] B ,l,i(p B r + p B .r))tp ( ' ±) (r) (16) 
For instance, if one considers i^+)(r), the quantity i?i(q) is given by 



and one has to evaluate 



H 1 (d) = (h,I 2 ,h,h) T 



(17) 
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1 f e~ r 

h = —j= / drcxp(iq.r) xF^irjE, 1, i(p B r + p B .r)) 

v 71 " J r 

To do that, we introduce the well-known integral [16] 



J(A) = / rfrexp(zq.r)— — iFi(w? s , 1, i(p B r + p B .r)) 
The other integrals can be obtained by noting that 



47T 



<? 2 + A 2 



exp 



ir\B ln( 



q 2 + X 2 



q 2 + A 2 + 2q.p B - 2iAp s 



and 



cos^exp(iq.r) = — — — cxp(zq.r) 

r oq z 



sin#cxp(z</>) exp(iq.r) = — -(tt — h exp(iq.r) 

r dq x dq y 



The second term in the transition amplitude given in Eq. (13) is 



c(2) _ o(2),l , < ? (2),2 



(18) 



(19) 



(20) 



(21) 



(22) 



with 



c(2),l 
6 /i 



/ 



u(p f ,s f ) u(j)i,Si) 1 u(p Bl s B ) Q E B , , r±), . .. .. 
dr ,^ 7^-^== __ _ 7 2) [ 7 2) _ -^]^(r)i^i(«7B + l,2,t(p B r + p B .r)) 



8tt 2 



x exp(-zp B .r)[exp(iA.r) - 1] exp(7T77 B /2)r(l - ir/ B )—5(E f + E B - E, t - e b ) 



and 



5(2), 2 



u(p f ,Sf) o u(pi,Si) 1 u{p B ,s B ) , (±) 



VW^ 7 w 2c V2SbF 

p_B-r)[exp(zA.i 

In Eq. (24), <p'(+)(r) for spin- up is given by 



87T 

x exp(-ip B .r)[exp(zA.r) - 1] cxp(7r?7 B /2)r(l - ir] B )—S(E f + E B - E t - e b ) 



(23) 



(24) 



</+)(!•) = N L 



( i/2c \ 


cos(<9) 
V sin(6»)e^ J 



Using the standard procedures of QED [19], one obtains for the unpolarized TDCS 



da 



9 X/ o X/ 



dE B dfl B dfl f 2 ^ ^ 2 ^ dE B dfl B dfl f 

1 Si.Sf S B S t 1 



(25) 



(26) 



evaluated for = Ei + e& — £3, where J2 St denotes the averaged sum over the spin states of the target atomic 

hydrogen with 



da 



1 Nbgj CXp(7T77 B ) 2 
^7 — — |1 (1 - IT}B)\ 



dE B dn B dfl f 64c 6 tt 3 \ Pi \ A 4 



(27) 



To our knowledge, in the expressions of S^'' and S,/' , a new integral occurs. We have calculated this integral 
analytically. Details of its derivation are given in the Appendix. This integral is 

J(A) = J drcxp(zq.r)^ t F 1 (i Vb + 1, 2, i{p B r + p B .r)) = 2*1 (m7b + 1, 1, 2, -2 (q ' P f 7 ) (28) 



(g 2 + A 2 )- 



g 2 + A 2 
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All the calculations in Eq. (27) can be done analytically and only five terms from nine are non zero, the diagonal 



terms 



5(1) 
5 /» 



S 



(2),1 
fi 



s 



(2) ,2 



and S^S^'' 1 as well as Sy/^'SY/ . In Eq. (26), the different sums over spin states 



5(2),ltc(l) 



give rise to the following results 

2 Ss 



«(P/>s/)7(°i)«(Pi>^) =2c 2 ( 

Bi 

(2) c 



2EiE, 



(Pi-Pf) 



|E st (--) = i 
r 



(29) 



III. RESULTS AND DISCUSSION 



A. Coplanar asymmetric geometries 

We begin our discussion by considering well known re- 
sults in the non relativistic domain, namely the results 
of Byron and Joachain [17] and those of Berakdar [21]. 
All these results are obtained in the coplanar asymmetric 



Experimental date 
SRC E3 A 
N RGB A 




FIG. 1: The two TDCSs. The solid line represents the rela- 
tivistic TDCS in the semi-relativistic Coulomb Born approx- 
imation, the long-dashed line represents the corresponding 
TDCS in the non relativistic Coulomb Born approximation. 
The incident electron kinetic energy is T = 250 eV and the 
ejected electron kinetic energy is Tb = 5 eV . Experimental 
data is from 1241 



geometry. Let us consider the process whereby an inci- 
dent electron with a kinetic energy Ei = 250 eV scatters 
with a hydrogen atom. The ejected electron is observed 
to have a kinetic energy E B = 5 eV and the scattered 
electron is observed having an angle 0/ = 3. In this 
particular case, the CBA is not as accurate as the results 
obtained within the framework of the Eikonal Born series 
[17] which contains higher order corrections. 

Nevertheless, as it can be seen in Fig. 1, the agreement 
between the non relativistic and semi-relativistic results 
is good since we obtain two identical curves. However, 
even in this non relativistic regime, small effects due to 




-1 SO -90 



Angle e B (degree) 



FIG. 2: The ratio T DC S{SRC B A) /T DC S{N RC B A) as 
function of the angle 9b with Of — 3°. The incident elec- 
tron kinetic energy is T = 250 eV and the ejected electron 
kinetic energy is Tb — 5 eV 



the semi-relativistic treatment of the wave functions we 
have used, show that there are indeed small effects that 
can only be tracked back to the spin. Indeed, if we plot 
the ratio of the semi-relativistic TDCS and the non rela- 
tivistic TDCS, it emerges that however small, these spin 
effects can reach 0,45% for some specific angles. We re- 
call that the TDCS has extrema, in particular when the 
direction of pe coincides with that of the vectors A and 
—A and this can be seen in Fig. 2. 

In the former case, the extremum is always a maximum 
and in the latter case the extremum is a local maximum. 
The two TDCSs exhibit in this geometry a forward or 
binary peak with a maximum in the direction of A and 
a recoil or backward peak in the opposite direction — A. 
The locations of such extrema are 9 B w —128° with a 
ratio equal to 1.00234 and 9 B w 52° with a ratio equal 
to 1.00185. These mechanisms for the emergence of the 
binary-recoil peak structure are also present even when 
one uses the simplest description in which plane waves 
for incoming and outgoing particles are assumed [22]. 

Now, if we compare our result with the result obtained 
by Berakdar [21], we also obtain a good agreement. But 
before beginning the discussion proper, let us recall the 
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formalism used by Berakdar. His calculations were per- 
formed within a model where the three-body final state is 
described by a product of three symmetrical, Coulomb- 
type functions. Each of these functions describes the mo- 
tion of a particular two-body subsystem in the presence 
of a third charged particle. Thereafter, he made a com- 
prehensive comparison with available experimental data 
and with other theoretical models. He ended his study by 
concluding that generally, good agreement is found with 
the absolute measurements but that however, in some 
cases discrepancies between various theoretical predic- 
tions and experimental findings are obvious, which high- 
lights the need for a theoretical and experimental bench- 
mark study of these reactions. In Fig. 3, we compare 




Angle i|> B (degree) 




Angle t> B (degree) 



FIG. 4: The ratio TDCS(SRCBA) /TDCS(NRCBA) as 
function of the angle 9 b with Of = 357°. The incident elec- 
tron kinetic energy is Tj = 250 eV and the ejected electron 
kinetic energy is Tb = 5 eV 



B. Binary coplanar geometries 

The relativistic regime can be defined as follows : 
when the value of the relativistic parameter 7 = (1 — 
(/3/c) 2 ) -1 / 2 is greater that 1.0053, there begins to be 
a difference between the non relativistic kinetic energy 
and the relativistic kinetic energy. This numerical value 



FIG. 3: The two TDCSs. The solid line represents the rela- 
tivistic TDCS in the semi-relativistic Coulomb Born approx- 
imation, the long-dashed line represents the corresponding 
TDCS in the non relativistic Coulomb Born approximation, 
the symbols square and circle respectively represent the for- 
malism of DS3C and the experimental data. We keep the 
same energies as in Fig. 1. Experimental data is from |f| 



our results with those obtained by Berakdar for an inci- 
dent electron kinetic energy = 250 eV for the case of 
a coplanar asymmetric geometric where Of = Bb = 90°. 
The ejected electron kinetic energy is Eb = 5 eV and 
4>f = 357°. What is remarkable is the agreement be- 
tween our results and his bearing in mind that he used 
the DS3C formalism (DS3C stands for dynamical screen- 
ing theory with three Coulomb- type functions). Another 
atypical result related to our calculations is the behavior 
of the ratio of the TDCS{SRCBA)/TDCS(NRCBA) 
where now the maxima of this ratio correspond nearly to 
the local minima of the TDCS when plotted as a function 
of the angle 4>b- This is shown in Fig 4. However, there 
is no rule that can be inferred from the behavior of this 
ratio since when performing various simulations even in 
the coplanar asymmetric geometry but with increasing 
values of the incident electron kinetic energy, there are 
many regions not close to the binary or secondary peaks 
that present maxima or minima. 



A 

A 


SRCBA 

NRCBA 

RPWBA 









Angle Q B (degree) 

FIG. 5: The three TDCSs scaled in 10" 3 . The solid line rep- 
resents the relativistic TDCS in the semi-relativistic Coulomb 
Born approximation, the long-dashed line represents the cor- 
responding TDCS in the non relativistic Coulomb Born ap- 
proximation. The short-dashed line represents the rlativistic 
plane wave Born approximation. The incident electron kinetic 
energy is Ti — 2700 eV and the ejected electron kinetic en- 
ergy is T B = 1349.5 eV and 0/ = 45° 

of the aforementioned relativistic parameter corresponds 
to an incident electron kinetic energy of Ei — 2700 eV. 
Because there is no experimental data available for this 
regime, we simply compare our results with those we have 
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previously found when we introduced the RPWBA |23( 
(relativistic plane wave Born approximation)to study the 
ionization of atomic hydrogen by electron impact in the 
binary geometry. In Fig. 5, it is clearly visible that the 
three models (NRCBA, SRCBA and RPWBA) give the 
same results which was to be expected since in this geom- 
etry, the use of a Coulomb wave function is not necessary. 



- SRCBA 
- NRCBA 




Angle e B (degree) 



FIG. 6: The two TDCSs scaled in 10" 5 . The solid line repre- 
sents the relativistic TDCS in the semi-relativistic Coulomb 
Born approximation, the long-dashed line represents the cor- 
responding TDCS in the non relativistic Coulomb Born ap- 
proximation. The incident electron kinetic energy is Ti = 
25000 eV and the ejected electron kinetic energy is Tb = 
12499.5 eV and f =45° 



In Fig. 6, there is a shift of the maximum of the TDCS 
in the SRCBA towards smaller values than Ob = 45° and 
this remains the case for increasing values of the kinetic 
energy of the incident electron. The origin of this shift 
stems from the fact that the main contribution to the 
TDCS comes from the term Hi(q) given by Eq. (16). 
This term contains a dominant integral I\. When plot- 
ting the behavior of I\ as function of the angle 9b , and 
with increasing values of Ei, one observes the shift we 
have mentioned as well as the fact that in the relativis- 
tic regime, the TDCS(SRCBA) is always lower than the 
TDCS(NRCBA). 



IV. CONCLUSION 



In this work, we have developed an exact semi- 
relativistic model in the hrst Born approximation that is 
valid for a wide range of geometries, simple in its mathe- 
matical structure and that allows to find previous results 
using sophisticated non-relativistic models. This model 
gives good results if the condition Za <C 1 is fulfilled. 



APPENDIX A: ANALYTICAL CALCULATION 
OF THE INTEGRAL J(A) 

Before turning to the analytical calculation of the inte- 
gral J(A) proper, let us recall how the integral /(A) [16] 
can be obtained. This is explained without any detail 
in [20]. Using parabolic coordinates, one has to evaluate 
the following integral 



,-Ar 



= J dr exp(zQ.r) cxp(— ips.r) 

The choice of the scalar product Q.r chosen is [6] 
Q.r = ^Q(£ — T)) cos 7 — 4-\/^cos(^sin7 

Performing the various integrals, one finds 
2tt 



(Al) 



(A2) 



/(A) = 



A - i(Qcos7 +pb) 
x iFi(ii]B,l,ipB£) 



POO 

/ eZ£exp(-/0 
Jo 



(A3) 



We use the well known result [20] 

/ dtexp(-\t)iFi(a,l,kt) = A a-1 (A- k)~ a (A4) 
Jo 



with 
A = /! 



Q 2 sin 2 7 



2[X -i(Q cos -f + p B )} 2 



[A + i(Qcos7 + Pb)] 



and a = ir\B and k = ips- This gives the result : 

47T 

/(A) = (Q 2 + X 2 + P 2 B + 2Qp B cos 1 ) 

Q 2 + A 2 + p% + 2Qp B COS7, 



(A5) 



x exp 



irjB ln( 



A2 - Pb - 2i ^PB 



(A6) 



To recover the integral /(A) given in EQ. (19) of the text, 
one has to make the following substitutions : 



Q = q + Ps 

Qp B cos j = Q.p B = q.p B -P% 

It is then straightforward to find that 



(A7) 



Q 2 + A 2 +p% + 2Qp B COS7 = q 2 + A 2 
Q 2 + A 2 - p% - 2i\p B = q 2 + A 2 + 2q.p B - 2i\p B 

(A8) 

so that 



1(A) 



-in 



{q 2 + A 2 ) 
To calculate 
J(A) = 



exp ir] B M — — ) 

q z + A 2 + 2q.p B - 2i\p B 

(A9) 



dr exp(iQ.r) exp(— ip^.r)- 



,-Ar 



x iF\{irjB + 1,2, i{pB r + Ps-i")) (A10) 
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one uses the same procedures to obtain 
2tt 



J(A) = 



/ 

Jo 



A - z(Qcos7 +ps) 

X 1^1(^5 + 1,2,^0 

2tt 1 

A - z(Qcos7 +Pb) M 

x 2 F 1 {ii 1 B + 1,1,2,—) 
I- 1 



rf^exp(-^) 



Performing the various substitutions, one gets the follow- 
ing new (as far as we know) analytical integral 



(All) 



-Xr 



J(X) = / (irexp(iq.r) \Fi(irj B + 1, 2, i(p B r + p B .r)) 
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(<Z 2 + A 2 ) 



2F1 ir/s + 1,1,2,-2 



g 2 + A 2 



(A12) 



r 



We have tested this analytical result by performing the 
integral using two gaussian quadratures because we have 
assumed without loss of generality both q and p B to be 
parrallel to the Oz axis. The first one, a Laguerre gaus- 
sian quadrature (32 points) to integrate over the radial 
variable r and the second one, using a Legendre gaussian 
quadrature (32 points) to integrate over the angular vari- 
able 9. The agreement between the analytical result and 
the numerical result is excellent. To illustrate this point, 
we give as an exemple the results obtained by the two 



methods for the following random values of the relevant 
parameters. For A = 1 , | q| = 1 .015055, |p| =0.1055098. 
The exact result is 

JexactW = (0.57355896,0.12458510) (A13) 

and the numerical result is : 

Jnum(A) = (0.57355899,0.12458507) (A14) 
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